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(Nl Abstract 

We present an elementary derivation of the width of photon decay in a weak mag- 
netic field using the semiclassical method of worldline instantons. The calculation is 
O | generalized to a model of quantum electrodynamics with broken Lorentz symmetry. 

Implications for the search of deviations from Lorentz invariance in the cosmic ray 
I, experiments are discussed. 

t> 1 Introduction 
o 
o 
m 



Semiclassical methods are widely used in modern quantum field theory. They provide 
a powerful tool to investigate non-perturbative phenomena. The well-known example is 
the false vacuum decay [1, 2]. Equations of motion have a solution, called bounce, that 
interpolate between false and true vacua. The probability of decay is proportional to the 
exponent (with the minus sign) of the Euclidean action evaluated on the bounce. 

The similar calculation arises in another class of processes, particle production in ex- 
ternal backgrounds. The simplest example is the Schwinger effect — spontaneous creation 
of electron-positron pairs in external electric field. The probability of pair production is 
expressed through the thermal partition function of a certain quantum mechanical problem 
[3]. This partition function is evaluated in the saddle point approximation; the solutions of 
the saddle point equations, called 'worldline instantons', are interpreted as trajectories of a 
particle in an auxiliary periodic time. Semiclassical treatment of Schwinger-like processes 
has been generalized to the cases of time-dependent and inhomogeneous electromagnetic field 
[4] and to the photon-stimulated Schwinger pair creation in application to the laser physics 
[5]. In this article we study the similar process — the photon decay into an electron-positron 
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pair in magnetic field. We will be interested in the weak-field limit where the semiclassi- 
cal method is applicable. For a review of methods used in the opposite case of the strong 
magnetic field see [6]. 

The probability of the photon decay in magnetic field was calculated long time ago 
independently by [7, 8]. The corresponding matrix element has been computed in the semi- 
classical approximation in terms of the overlap of electron wavefunctions (eigenfunctions of 
the Dirac equation in the uniform magnetic field) in the coordinate representation. In the 
weak-field limit the photon decay is exponentially suppressed, 



Here oj is the energy of the photon, H is the value of the uniform magnetic field, (p denotes 
the angle between the photon momentum and the magnetic field, e and m are the electron 
charge and mass; it is assumed that the photon energy is much higher than the electron 
mass, usimp ^> m (but still ueHsiwp <C m 3 ). The calculation of [7, 8] is technically quite 
involved. The approach based on the 'worldline instantons' adopted in the present paper is 
significantly simpler and has a clear geometrical interpretation. To the best of our knowledge, 
it has not been applied to the process of pair production in magnetic field so far. 

Due to its simplicity, our method can be easily generalized to models beyond the standard 
QED. This is illustrated in the second part of the paper where we study the process of gamma 
decay in magnetic field in a model of electrodynamics without Lorentz invariance (LI). It is 
possible that deviation from LI can appear at very high energies, unaccessible to present-day 
accelerators. This scenario is suggested by several approaches to the theory of quantum 
gravity [9, 10, 11, 12, 13] and the energy scale where deviations from LI become significant 
is naturally assumed to lie at the Plank mass M p or a few orders below. 

Remarkably, this type of Lorentz violation (LV) can be constrained by cosmic ray obser- 
vations, see [14] for recent review. Indeed, energies attained by particles of ultra-high energy 
cosmic rays (UHECR) greatly exceed those achieved in laboratory. Reaching the Earth, 
cosmic ray primaries interact in the atmosphere and create showers of descendant particles, 
that can be detected experimentally. The characteristics of the shower depend on the type 
and energy of the primary. 

The process of photon decay in magnetic field plays an important role in these con- 
sideration. When an UHECR photon (with energy above 10 19 ' 5 eV) reaches the magnetic 
field of the Earth, it decays into an electron and positron. The latter, in turn, produce 
photons by the synchrotron radiation, giving rise to an electromagnetic cascade. In this 
way the magnetosphere shower called preshower [15, 16] is created at the altitude of several 
hundred kilometers above the Earth surface. Preshower accelerates the subsequent shower 
development in the atmosphere and provides a unique signature for photon-induced show- 
ers: pair production probability depends on the perpendicular component of the magnetic 
field and hence on the arrival direction of the photon. We will find that possible LV signifi- 
cantly affects the probability of the preshower formation. Thus, experimental detection of a 
preshower would provide a sensitive probe of LV. 




(1) 
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The paper is organized as follows. In Sec. 2 we derive Eq. (1) by the semiclassical method 
of worldline instantons. In Sec. 3 we calculate the width of photon decay in magnetic field 
in QED with LV. Sec. 4 is devoted to discussion. 

2 Photon decay in a weak magnetic field in standard 
QED 

Consider a photon with four-momentum = (w,k) propagating in the uniform magnetic 
field H at an angle p> to its direction. We choose the coordinate system where the magnetic 
field points along the a>axis, H = (if, 0,0), and the spatial photon momentum lies in the 
(x, y)-plane, k = (u cos ip, u sirup, 0) . The photon decay into an e + e" pair is kinematically 
allowed if u> sin p > 2m. We will assume the photon energy to be well above this threshold, 
usinp 3> m. 

To find the rate of the photon decay we adopt the method similar to that used in [3, 4, 5] 
for the semiclassical analysis of the Schwinger process. It was shown in these works that 
in the leading approximation the answer is insensitive to the spin of the electron. As in 
the present paper we are interested only in the leading order result, we choose to work for 
simplicity with the scalar QED, described by the Lagrangian 1 

C = -\f, u F^ + D^D^ - m 2 0*0, (2) 

where the covariant derivative is defined in the usual way, D^ifi = (d^ — ieA^) <fi. 

It follows from the optical theorem that the rate of photon decay is proportional to the 
imaginary part of the polarization operator 

T = ^{k)e v {k)\mll^{k), (3) 

where e^(k) is the photon polarization vector which we choose to be real. As usual, 

is given by the Fourier transform of the correlator of two electromagnetic currents j M = 

n, v (k) = j d A y ^{j,{y/2)U- y /2)) . (4) 

In the leading order one can neglect the contribution of virtual photons into this correlator 
which is then expressed in terms of the partition function of the charged scalar in external 
electromagnetic field, 

0,fa/2)>(- y /2)) = z^ ^naA.L/*) ™ ' <5) 

1 We take the signature (+, — , — , — ) for the Minkowski metric. 
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where 



Z[A,} = J D<l>*D<t>e-S d4 <\ D ^ 2+m2 ^) =det (-Dl + m 2 ) = expTrln (-Dj + m 2 ) . (6) 

Note that in defining the partition function we have performed the Wick rotation to the 
Euclidean signature. 

At the next step we use the formula: 



-InM 



[ % (e- MT ~ e- T ) 



This leads to the expression for the partition function in terms of the integral over the 
"proper time" T. 

r x , ; ' / ■ 

-m?Tr 



Z[AJ[ = Z exp 



00 dT 



T 



Tr 



(e™S) 



(7) 



The operator (— Dfy can be interpreted as the quantum-mechanical Hamiltonian of a point 
particle in four-dimensional space. Thus, Tr can be considered as its thermal 

partition function with the proper time T playing the role of the inverse temperature. It is 
more convenient to work in the Lagrangian formalism, so we make the Legendre transfor- 
mation and consider the functional integral representation: 



Tr 



(e T ^) = 



Dx^e 



- Io dT (-r+ ie± ^ A ^j 



p.b.c. 



Here we have introduced an auxiliary time r and the notation p.b.c. means periodical bound- 
ary conditions x^O) = x^(T). The partition function becomes 



Z[A^ = Z exp 



dT 



-m z T 



b.c 



- fo dT (^r+ ie± ^ A ^) 



(8) 



We now return to the polarization operator (4). Each variational derivative of the parti- 
tion function with respect to produces an insertion of the combination § drx^{T)5{x{T) — 
y) in the functional integral. Also we rescale the auxiliary time r to make it vary from to 
1. This gives, 

Un(y/2)M-y/2)) oc f ^ / Dx^ I dn I dT 2 x lx {T 1 )x 1/ {T 2 ) 

JO 1 Jp.b.c. J J 



The width of the photon decay is given by the imaginary part of the polarization opera- 
tor in the momentum representation ImU^^k) = Im J d 4 y e^U^ly). Taking the Fourier 
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transform for the imaginary part of the polarization operator 2 , one obtains 



where 



ImU^k) oc Im f D Xfl I dr x I dr 2 ^(ri)x^(r 2 )e S ™^ T ^\ 

JO 1 Jp.b.c. J J 

S m [x^ n, r 2 ] = m 2 T + J dr + ieA^x^j - ik^ (x^n) - a^(r 2 )) . 



(9) 



(10) 



This expression has the form of the Euclidean action of a relativistic particle in the external 
electromagnetic field with two sources of opposite signs located at the proper times T\ and 
t 2 . The strength of these sources is determined by the photon momentum. Integrating out 
the parameter T one can obtain the standard form of the relativistic particle action. 

We will evaluate the r.h.s. of (9) in the saddle-point approximation. To this end, firstly 
we find the saddle equations for T and x M (r). Their solution gives the saddle-point classical 
trajectory x^(r). At the second step the trajectory is substituted into the action (10). Let 
us fix the gauge = —^F^ u x u . Varying over x^, we obtain (separately for the time and 
space components of x^), 



x 
2T 

Xi 

2T 



= u (5 (r - ti) - 5 (r - r 2 )) , 

— ieFijij = —iu [5n cos (p + 5 i2 sin ip] (5 (r — Ti) — 5 (r — r 2 )) 
The variation over T yields, 



m 2 — 



AT 2 



0. 



(11) 
(12) 

(13) 



We are looking for a solution of Eqs. (11)— (13) that describes a closed trajectory in 4- 
dimensional spacetime. Solutions of this type are called "worldline instantons". Note that 
in general they can be complex (cf. [5]). The solution exists if \t± — r 2 | = \. Without loss 
of generality we set t± — 0, r 2 = |. Our solution is composed of two hyperbolic arcs (see 
Fig. 1) defined on the segments r G (0, \) and r G (|, 1) respectively, 



for < t < 



x, 



cl 



ojT[t- 



71 

xf = — %Ar\ cos ip ch.— 



for - < t < 1 : 



x% = -uT [ T 



2>- 



xf = iArj cos p ch- 



x 2 = iA sin ip sh I rj ( r 



xf = A sin (p 



ch ( 77 ( r 



x 2 = —iA sin p sh ^77 — ^ ^ , xf = —A sin p ch ^rj — ^ 



— ch- 
4 

(15) 



ch- 



2 Strictly speaking, the integral here should be performed over the configurations satisfying x(t\) = — .t(t 2 ). 
We will ignore this restriction because it does not affect the final result. 
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Figure 1: 'Worldline instanton' — the classical trajectory describing pair production in the 
magnetic field for the case when the photon momentum is perpendicular to the field direction, 
ip = tt/2. The projection of the trajectory on the plane (1x2, £3) is shown. 



Here parameters A, r] are determined from Eqs. (11), (12): 

CO 

Substituting the solution (14)-(16) into (13) one obtains 3 : 

m 4m 



(17) 



ueH sin (p 



The next step is to evaluate the action (10) on the solution. After a straightforward 
calculations we obtain: 

S [ x <*] = g f . (18) 

The semiclassical method is valid as long as the classical action is large 4 , S[xr[] ^> 1. Com- 
bining everything together, we obtain the photon decay width 5 : 

where is the derivative of the classical solution 6 and J\f is a pre-exponential factor com- 
ing from the integration over fluctuations near the classical solution (14)-(16). Using the 



3 This formula is valid in the regime wsin</> ^> 2m. The exact expression reads T = ^ arcth UJ 2 s ^ 1(p - 

4 As an example, let us consider the geomagnetic field H ~ 0.5 G and take sin<^ ~ 1. Then the method is 

applicable for photons with the energy to < 10 19 5 eV. 

5 In general, one should sum over contributions of all possible classical solutions, we consider only the 

dominant one. 

6 Note that the discontinuity of x„ at t = 0,1/2, implied by Eqs. (11), (12), is proportional to and 
thus vanishes when contracted with the polarization vector. 
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reasoning similar to [3] one can show that the prefactor J\f has a single negative mode, corre- 
sponding to the variation of the size of the worldline instanton. Therefore, according to the 
standard arguments [2] the prefactor is imaginary and the decay width is nonzero. In our 
approximation we are interested only in the leading exponential behavior. Thus we neglect 
the prefactor and arrive at the formula (1). 



3 Generalization to QED with Lorentz violation 

The simplicity of the method presented in the previous section allows us to easily generalize it 
to non-standard theories. In this section we analyze the sensitivity of the photon decay width 
in the magnetic field to possible deviations from Lorentz invariance. To make the calculation 
concrete we need to specify the model. We consider the analog of the model introduced in 
[17], where we replace the spinor electron field by a charged scalar. The Lagrangian reads, 

C = - l -F^ v + D^*D^ - m 2 0*0 + (19) 

+ -F ki (-2x+ F k > - ^-D 2 6*DU. 
4 k] V M 2 J M 2 

Here the first line represents the standard LI scalar QED, while the second line contains 
extra LV operators of dimension 4 and 6; x, £ and g are dimensionless coefficients, M is a 
parameter of order the Planck mass. The form of the Lagrangian is fixed by requiring the 
theory to be rotationally invariant in the preferred frame, gauge invariant, and CPT- and 
P-even. The full list of restrictions on the theory and their motivation are discussed in [17]. 

From (19) one obtains the dispersion relations for photons and (scalar) electrons and 
positrons, 

7 : c 2 = A; 2 (l-2x) + |^, (20) 

e ± : E2 = ™ 2 +P 2 + 2 j^- (21) 

These differ from the standard case, and as a consequence the kinematics of various reactions 
is modified. In particular, the reactions relevant for propagation and detection of the cosmic 
rays are affected which makes the cosmic ray experiments sensitive to LV. 

As discussed in [17], in general there are other consequences of LV that may also be 
important. Thus, the sums over polarizations entering the calculation of reaction rates 
are also modified; new interaction vertices appear from the last term in the Lagrangian 
(19). However, these modifications affect only the pre-exponential factor in the polarization 
operator. As we are interested only in the leading exponential behavior, we neglect this type 
of correction in what follows. 

The optical theorem is based only on unitarity and does not rely on LI. Thus, to calculate 
the photon decay rate in magnetic field we can still use the formula (2). As before, we neglect 
the contributions of virtual photons, so that the polarization operator Yl^ v is again related 
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to the second variational derivative of the partition function Z[A^\. On the other hand, the 
latter is modified due to LV. In the proper time representation it reads, 



Z\Ay\ = Z exp 



dT 



e m2T Tre T < 



Di 



. 2 3 
M 2 



(D 2 ) 2 



(22) 



where we have rotated to the Euclidean time. Here the additional term — ^(D 2 ) 2 in the 
inner exponent comes from the four-derivative electron kinetic term in the second line of (19). 
The next steps are the same as in Sec. 2. One interprets the trace in (22) as a quantum 
mechanical statistical sum and writes the functional integral representation for it in terms 
of the point-particle action. This leads to the expression (9) for the imaginary part of the 
polarization operator, where now 



3,i 



m 2 T + 



dr 



AT ie ^ Xtl + 8T 3 M 2 j 



ik^ (av(ri) 



x, 



(r 2 )) • 



(23) 



In deriving this formula we have assumed, 



gxt 



M 2 T 2 



< 1. 



(24) 



We will see later, that this condition is equivalent to the requirement that the LV corrections 
to the electron dispersion relation (the last term in (21)) is small compared to p 2 . Note that 
this still allows the LV correction to be of order or larger than the (squared) electron mass 



m . 



I ■ 2 \2 

Eq. (23) differs in two respects from the LI case (10). First, the term ^fgjp appears in 
the point-particle Lagrangian. Second, the momentum k^ of the initial photon must satisfy 
dispersion relation (20). 

We now evaluate the integrals over T and over x M in (9) by the saddle point method. 
Varying the action (23) over T we obtain, 



m 



2 "V 



3g(x 



2\2 



4T 2 8M 2 T 4 ' 



(25) 



Variation over x M gives the equations of motion. The time- component of the equations does 
not change compared to the LI case, see Eq. (11). On the other hand, the spatial equations 
get modified, 7 



Xi 

2T 



1 + 



M 2 T 2 



-ik [8n cosy? + 5 i2 smtp\ (S(t - Ti) - S(r - r 2 )) . (26) 



The solution of Eqs. (11), (26) has the form (14), (15), but with different parameters A and 
V, 



A = 



k 



2eHch\ ' 



r] = 2TeH 



1 + 



7 Note that on the l.h.s. of (26) we have omitted the term — — 



M 2 

2T 3 M 2 -frf , as x'j is the integral of motion. 



(27) 
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Substituting this into eq. (25) and solving it with respect to T we obtain, 



eH sin tp 

Note that none of the terms under the square root can be neglected. Coming back to the 
condition (24), it is now straightforward to check that on the solution it reduces to 

- — < 1. 

M 2 

As advocated before, this is nothing but the requirement that the ratio between the third 
and the second term on the r.h.s. of (21) is small for the real electrons (positrons) produced 
in the photon decay. Clearly, this is satisfied for all astrophysically relevant photon energies 
uj if g is not much bigger than 1. 

Next we substitute the solution (14), (15), (27), (28) into the action (10) and obtain, 

/ 2 \ 3 / 2 

uj 2 If 2m\ goo 2 £co 2 



s M = m^UJ + fKp-w + 2 "\ ■ < 29 > 

Following [17] we introduce the combination uj^v that characterizes the LV contribution into 
the kinematics of the reaction, 

£cu 3 guj 3 . , 

^ = -xc+ — -— . (30) 

Using this notation the width of the photon decay is cast into the form, 

8m 3 / uj • ojlv\ 3 / 2 



T oc exp 



3ujeH sin ip 



/ _ uj ■ ujlv Y 
V 2m 2 J 



This is the main result of this section. Let us analyze it. Even small negative ujlv ^ — 



exponentially suppresses the width of the photon decay in magnetic field. On the other 
hand, even small positive ujlv ~ ^j- decreases the absolute value of the exponent and the 
decay becomes unsuppressed (of course, the semiclassical approximation breaks down in this 
case, cf. the discussion at the end of Sec. 2). 

Our result admits the following interpretation. Consider for simplicity the case when 
LV is present only in the electron sector, x = £ = 0. Introduce the effective momentum- 
dependent electron mass by the formula 

m 2 eff (p)^E 2 - P 2 = m 2 + 2 -l^. (31) 

In terms of this notation the formula for the pair-production width takes the standard form 
(1) with m replaced by m e ff(uj/2) — the effective mass of the produced electron (positron). 
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The larger the effective mass, the more suppressed is the photon decay, and vice versa. One 
concludes that in the leading approximation the effect of LV on this process is completely 
encompassed by the kinematics. 

Before finishing this section let us discuss how our calculation is modified in the case 
of a more general pattern of LV. Generalization to an arbitrary LV in the photon sector is 
straightforward: the only change amounts to a different relation between u and k in Eqs. (11) 
and (26). The form of the equations remains the same, and by literally repeating the above 
calculation one can analytically find the suppression exponent in a model with arbitrary 
dispersion relation u{k). 

More technical difficulties arise if we allow for an arbitrary dispersion relation for electrons 
(and the same for positrons). The saddle point equation (25) changes and, in general, cannot 
be solved analytically: all terms in this equation are comparable, so we cannot perform an 
expansion in small parameter. Thus, for a general electron dispersion relation one can find the 
width of photon decay only numerically. However, it appears on the physical grounds that, 
at least qualitatively, the leading effect of LV will be again encompassed by the kinematics. 
Thus, a faithful estimate of the suppression exponent can be obtained by substituting the 
effective mass m e ff defined by the first equality in (31) in the standard formula (1). 

4 Discussion 

We have studied the process of photon decay into an electron-positron pair in external 
magnetic field by the semiclassical method of worldline instantons. A technically simple 
derivation of the leading exponential behavior of the width of this process has been presented. 

We have shown that the method can be easily generalized to the extension of QED 
including possible deviations from the Lorentz invariance and illustrated this by an explicit 
analytical calculation in the model of LV scalar QED with dispersion relations quartic in 
momentum. We also discussed how this calculation can be in principle generalized to the 
case of arbitrary dispersion relations. It was found that the width of photon decay in weak 
magnetic field is exponentially sensitive to the LV contributions. 

Let us discuss the implications of our results for the tests of LV involving cosmic ray 
observations. To date, no UHE photons with energy 10 19 ' 5 eV or more have been detected. 
However, there are reasons to expect non-zero flux of photons at such high energies. The 
break in the spectrum of cosmic ray hadrons at energy ~ 10 19 ' 6 eV has been detected in- 
dependently in three experiments [18, 19, 20]. If the dominant fraction of the cosmic ray 
primaries are protons this break is naturally identified with the GZK cut-off [21, 22]: sup- 
pression of the proton flux due to the interactions with the cosmic microwave background. 
In this interaction multiple charged and neutral pions are produced. Neutral pions, in turn, 
decay to photons called 'cosmogenic', or 'GZK' photons [23, 24]. The cosmogenic photons 
may be detected within the current decade by the Pierre Auger Observatory [25]. 

Imagine now the situation that a non-zero flux of photons with energies u ~ 10 20 eV 
has been observed. Imagine moreover that the threshold for the preshower formation in the 
Earth magnetosphere has been measured to coincide with the predictions of the standard 
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Lorentz invariant QED. This would imply that the quantity ulv defined in (30) satisfies 
\ojlv\ < 2m 2 /oj, or numerically \ojlv\ < 10~ 8 eV. Barring accidental cancellations between 
various terms entering (30) and taking M = 10 19 GeV this would translate into the stringent 
bounds |x| < 10~ 28 ; [£[, \g\ < 10~ n . The ability to constrain £, g at the level well below 1 
means that the preshower formation by UHECR photons is sensitive even to trans-Planckian 
breaking of LI. 
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